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We use high dimensional bosonization to derive an effective field theory that describes the Pomer- 
anchuck transition in two-dimensional Fermi liquids. The bosonization approach explicitly retains 
all low-energy degrees of freedom of the system. The resultant theory has dynamical exponent z — 2 
at tree level and upper critical dimension d c — 2, thus in 2D the system is at the upper critical 
dimension. These results differ from those of an earlier study based on integrating out fermions. 
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Introduction and Motivation — Quantum phase transi- 
tions [1] are of tremendous current interest to the physics 
community in general and condensed matter physicists in 
particular. Significant progress has been made in theoret- 
ical studies of quantum phase transitions in various spin 
and boson models [1]; examples include various magnetic 
transitions in insulating spin systems, and superfluid- 
insulator transitions in bosonic systems. Historically 
however, the first examples of quantum phase transitions 
studied theoretically are ferromagnetic and antiferromag- 
netic transitions in itinerant electron systems, or metals 
[2]. Ironically, quantum phase transitions in itinerant 
electron systems turn out to be among the most difficult 
to study theoretically, and our current understanding of 
these transitions is much more limited compared to those 
of spin and boson systems. The origin of the difficulty is 
not hard to see. Just like classical phase transitions, the- 
oretical studies of quantum phase transitions are based 
on (quantum versions) of Ginsburg-Landau type of free 
energy functional, and their analyzes based on Wilso- 
nian renormalization group (RG). The Ginsburg-Landau 
free energy is a functional of the order parameter, which 
are bosonic degrees of freedom; fluctuations of their long- 
wave length, low-frequency/energy components deter- 
mine the critical properties of the transitions. In itin- 
erant electron systems however, the fundamental low- 
energy degrees of freedom are the gapless electronic states 
near the Fermi surface; these fermonic modes at fi- 
nite wavevectors are difficult to incorporate within the 
Ginsburg-Landau- Wilson (GLW) paradigm. In the ap- 
proach pioneered by Hertz [2] and extended by Millis 
[3] and others, one decouples the electron-electron inter- 
action using Hubbard-Stratanovish transformation in an 
appropriate channel, integrates out the fermionic degrees 
of freedom, and arrives at a GLW-like free energy func- 
tional that involves the Hubbard-Stratanovish auxiliary 
field only, which are bosonic degrees of freedom and in- 
terpreted as the order parameter. It has been realized 
recently, however, that the procedure of integrating out 
gapless fermions leads to singularities in the expansion of 
the resultant bosonic free energy functional in terms of 



the order parameter, or its gradients [4]. Such singular- 
ities may invalidate the standard classification and ana- 
lyzes (based on power-counting) of various terms of the 
GLW-like theory obtained this way, and profoundly affect 
the critical behavior of the transitions [4] . The presence 
of such singularities is a consequence of integrating out 
gapless degrees of freedom, and also reflects the fact that 
the GLW-like theory obtained this way only provides an 
incomplete description of the low-energy physics of the 
system. At present there is no consensus on the appro- 
priate approach to study quantum phase transitions in 
itinerant electron systems in general. 

Recently, we studied the ferromagnetic transition in 
one-dimensional (ID) itinerant electron systems, using a 
different approach [5]. The idea was quite simple: In ID 
one can use the very powerful machinery of bosonization, 
which allows one to obtain a bosonic description of the 
fermionic system, without losing any fermionic degrees 
of freedom; i.e., one can represent the fermionic degrees 
of freedom using bosonic degrees of freedom faithfully, 
via bosonization. Extensive previous studies have estab- 
lished that the paramagnetic phase in ID is described 
by a free boson theory, known as the Luttinger liquid 
theory in this context; this is a consequence of the fact 
that all possible bosonic interactions are irrelevant at the 
Luttinger liquid fixed point, and scale to zero in the low 
energy limit. We have shown [5], on the other hand, at 
the Gaussian critical point of the ferromagnetic transi- 
tion, interactions are relevant, and must be included in 
the analyzes of the critical behavior, as well as the or- 
dered phase. We have performed such analyzes using 
RG, and found the behavior of the system at the fer- 
romagnetic critical point is quite different from that of 
the Luttinger liquid, due to the presence of interaction; 
this "non-Luttinger liquid" behavior is in close corre- 
spondence with the experimentally observed non-Fermi 
liquid behavior near magnetic transitions in higher di- 
mensions. It is worth emphasizing that the theory we 
developed using bosonization does not suffer from the 
singularities encountered in the Hertz-Millis approach in 
high dimensional systems discussed above, since no low- 
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energy fermionic degrees of freedom are lost [6] . 

The bosonization machinery has been generalized to 
higher dimensional systems [7-10]. In general bosoniza- 
tion is not as powerful in high D as in ID. This is be- 
cause to bosonize the fermions, one needs to divide the 
Fermi surface (which is a continuous manifold in high 
D in contrast to discrete points in ID) into many small 
patches that are almost fiat; the scattering processes that 
bring the electron from one patch to another show up as 
non-linear operators within bosonization, rendering them 
difficult to treat. What the bosonization approach does 
treat well are the forward scattering processes which leave 
the electrons in the same patch that they start from, as 
they are represented by quadratic terms of the boson 
field; these are precisely the Fermi liquid interactions. 
For this reason bosonization gives a very good descrip- 
tion of the Fermi liquid phase [7-10]. 

In the present work we use high D bosonization to 
study the Pomeranchuck transition in 2D spinless Fermi 
liquids, a quantum phase transition at which the Fermi 
surface develops anisotropy spontaneously, which is of 
considerable current interest [11-13]. The bosonization 
approach is particularly useful here because the Pomer- 
anchuck transition is driven by Fermi liquid interactions. 
The strategy here is similar to Ref. [5] . Our central result 
is the GLW theory for the transition, Eq. (16). From the 
action (16) we can readily read out that the tree-level 
dynamical exponent z = 2 at the transition, and upper 
critical dimension d c = 4 — z = 2. Thus in 2D the system 
is at the upper critical dimension. This differs from the 
result an earlier study [11] based on Hertz-Millis type of 
approach, which found z = 3 and d c = 1. We will discuss 
the origin of the difference of these two approaches. 

Bosonized Description of Fermi Liquids — We start 
by reviewing the bosonized description of Fermi liquids, 
which also sets the notation for the rest of the paper. 
As the first step we divide the Fermi surface into many 
small patches of linear size A <C kp, so that within each 
patch the Fermi surface is approximately flat, and for 
electronic states sufficiently close to the Fermi surface 
(within a cutoff distance A <C hp), the energy is approxi- 
mately linear in momentum. For each patch we introduce 
a patch density operator 

p(S,q) = ^0(S,k)4_ qCk , (1) 

k 

where S is the patch label, and #(S,k) is 1 when k is 
inside the box of size A D_1 A enclosing patch S (D is 
the dimensionality of the system), and zero otherwise; it 
ensures we are summing over states near patch S only. 
They satisfy the bosonic commutation relation 

[p(S, q),p(T, q')] - nas,T* q ,-q'(q ■ n s ), (2) 

with £1 = A D ~ 1 (L/2ir) D , and L is the linear size of the 
system and ns is the unit vector pointing in the outward 



normal direction of patch S. In the case of ID the Fermi 
surface patch reduces to a discrete Fermi point, and Eq. 
(2) reduces to the familiar ID commutator of density 
operators. 

The foundation of bosonization lies on the fact that 
both the kinetic energy and the Fermi liquid (or forward 
scattering) interaction terms of the electron Hamiltonian 
can be expressed in terms of p(S, q). Within the approx- 
imation that the kinetic energy is linear in momentum 
within each patch: es(k) w vp(S)ris ■ [k — kp(S)], where 
vp(S) is the Fermi velocity of patch S, the kinetic energy 
is quadratic in p(S, q): 

T = 5> F (S)/2n]p(S,q)p(S,-q). (3) 

S,q 

Non-linearity in the spectrum leads to higher order terms 
in p(S,q); they are irrelevant for the description of the 
Fermi liquid fixed point, but are important for the de- 
scription of the transition and will be discussed later on. 
The Fermi liquid (or forward-scattering) interaction takes 
the form 

^=^E y S,T(qMS,qMT,-q), (4) 

S,T,q 

where Vs.T(q) is the forward scattering matrix element 
between states in patches S and T; at the Fermi liquid 
fixed point its dependence on q is irrelevant and we can 
take the q = value in the long-wave length limit (here 
we do not consider possible singularity due to long-range 
interaction). The Fermi liquid Hamiltonian H = T + V , 
which is expressed exclusively in terms of the patch den- 
sity operators p(S,q), also has an equivalent Lagrangian 
description: 

L = L {p(S,ci)}-H{p(S,q)}, (5) 

where the dynamical term [9,13] 

L ° = E^( s < i)M s . -<0/fa ' «s) (6) 

S,q 

properly enforces the commutation relation (2). 

Two-dimensional Isotropic Fermi Liquid — We now 
turn our discussion to the case D=2, and assume isotropy 
so that the Fermi surface is a circle. In this case the patch 
S can be labeled by an angular variable 6, and in the 
limit A <C k F , 9 may be treated as a continuous variable 
ranging between and 2tt. In this case it is natural to 
perform a Fourier transformation with respect to 8: 

pM = ^Jd8e- lme p(9, q ); (7) 
P(0,q)= £ e me Pm (q). (8) 

m— — Trkp- /A 
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The cutoff in the range of to is due to the small but finite 
patch size A. We can now express Eq. (6) in terms of 
Pm(q): 

L Q = (ik F /2Afl) 

x J2 [3t/> m (q)]p-n(-q) I d9e« m -^ e /[ q -h(9)] 

q,m,n 

mk F ^ [dtp m (q)} , tie. 



2(-l)-5-e iW V-m-i(-q). (9) 



Afi ^ |q, 

q,m 11 odd J 

where # q is the angle of 2D vector q. It is very important 
to notice that in (9) p's with even to are coupled to p's 
with odd n only, and vice versa. This is because the 
integral over 9 vanishes when to — n is even, so in the 
last line of (9) the sum over I is for odd integers. 

On the other hand it is easy to show that both T and 
V are diagonal in to, and the Hamiltonian takes the form 
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m. , q 



(1 + F TO )p m (q)p_ m (-q), 



where the dimensionless Landau parameter 
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in which V(9) is defined through Vs.T(q) = V(9s—9t, q), 
and we take q = here in the long wave-length descrip- 
tion of Fermi liquid; the q dependence of V will become 
important and discussed later on. The stability of the 
isotropic Fermi liquid phase requires F m > —1 [9]. Thus 
combining (9) and (10) we obtain a quadratic theory in 
terms of complex bosonic variables p m (q) , that describes 
an isotropic 2D Fermi liquid. Notice that L involves one 
time derivative of p m (q); this indicates p TO (q) are con- 
strained degrees of freedom (or their conjugate momenta 
arc not independent variables). To bring L closer to the 
more familiar ID bosonized description of Luttinger liq- 
uids, we can integrate out p m (q) with odd m, and obtain 
the Lagrangian with even to degrees of freedom: 



T — T" 

J-'even — ^even 



(12) 
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even 
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m.n even 



e^-^+^dtP-U-qjdtp^d), (13) 



where A e = J^iodd 4(i+f,) and H even takes the same 
form of (10) but with sum over even m. Despite the fact 
that we have integrated out modes with odd m, L even is 
completely equivalent to L because L® ven contains two 
time derivatives and thus the conjugate momenta of p m 
are independent degrees of freedom, which actually cor- 
respond to p m with odd m. Alternatively we may choose 
to integrate out modes with even m to obtain a dual de- 
scription L dd] this is analogous to the case in ID where 



we may write down the Luttinger liquid Lagrangians in 
terms of either density or current fields, which are sym- 
metric and antisymmetric combinations of left and right 
moving fields respectively. Clearly L and L even are scale 
invariant field theories with dynamical exponent z = 1, 
because there is a factor l/|q| associated with each time 
derivative (or ui in frequency space) while # q is invariant 
under scale transformations. 

The Pomeranchuck Instability and Ginsburg- Landau- 
Wilson Theory for the Transition — Obviously the 
isotropic Fermi liquid phase becomes unstable when F m 
reaches and goes below —1 for any to; this is the Pomer- 
anchuck instability. Let us assume this occurs (without 
losing generality) in an even channel mo ^ 0. In this case 
we need to retain the next leading quadratic term in H 
for this channel, of the form 



V = a |q| 2 Pm (q)p-mo (-q) + 



(14) 



it originates from the q dependence of Vs ; T(q) m Eq. (4), 
neglected at the isotropic Fermi liquid fixed point due to 
its irrelevance there. Here we assume a > 0. To maintain 
stability, we also need to keep non-quadratic terms [13] 
of the form 

t' = b Y Pm(q)pn(q')p-m-n(-q - q') 

m,n, q,q' 

+c X! Pm (q)p»> (q')pj (q")p-m-n-«(-q - q' - q") 
+ ■■■, (15) 

whose origin is the nonlinearity of electron dispersion 
near the Fermi surface (b oc e"(kp) and c cx e'"(fcir)), 
again neglected at the isotropic Fermi liquid fixed point 
due to its irrelevance. To proceed we focus for the mo- 
ment on the mo channel where the instability occurs, and 
neglect its coupling to other channels. From Eqs. (13, 
14, 15) we obtain the following GLW effective field theory 
with Euclidean action 

S = J ^q^{(^l + |q| 2 +r)|^(q,c)| 2 

+ (-l) mo ^Re[e 2lm « e ^(q, w^-q, -co)]} 

+ uj d D xdr|0(x,T)| 4 + ---, (16) 

where 4> cx p mo is a complex bosonic field that plays the 
role of order parameter of the theory, the "mass" -like pa- 
rameter r cx 1 + F mo , and u cx c (also assumed positive) ; 
proper rescaling of the field as well as space-time coordi- 
nates have been performed to ensure the form taken by 
the quadratic terms in Eq. (16). This highly non-local 
action is somewhat similar to the one studied by Sachdcv 
and Senthil (see Eq. 4.11 of Ref. [14]) in a different con- 
text; the crucial difference here are the new quadratic 
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terms in the second line of Eq. (16), which encode the 
information about the symmetry properties of the order 
parameter as well as the Fermi surface dynamics. Clearly 
at the Gaussian critical point r = 0, we have dynami- 
cal exponent z = 2 from the u> and q dependence of the 
quadratic part of the action. The tree-level flow equation 
for the interaction is simply given by its dimensionality: 

fhi 

= (4- *-£>)«, (17) 

a logs 

where s the scaling parameter for spatial coordinates; we 
thus find the upper critical dimension of the theory to be 
d c = 4 — z = 2. Higher order couplings are irrelevant at 
D = 2. These are the central results of this paper. 

We now show that couplings to the other (non-critical) 
channels do not change the form of the action (16), or the 
associated critical behavior, for the following reasons, (i) 
For the other channels the action contains a non-zero 
mass-like term (1 + F m ) p m (q) p_ m (— q). Under the z = 2 
scaling the "mass" 1 + F m is a relevant operator that 
grows as s 2 under scaling [15]. Thus in the long-wave 
length limit all other channels become "infinitely mas- 
sive" and thus drop out of the low-energy effective action 
eventually, (ii) Before this limit is reached, one needs to 
integrate over modes in these channels with wave- vectors 
at the cutoff as RG proceeds; this generates various new 
couplings in the m channel, in addition to renormalizing 
the existing ones. Due to the presence of the mass-like 
terms, this is a well-behaved procedure, and due to the 
nature of the coupling in Eqs. (12,13) between the mo 
and other channels, the only possible singular new cou- 
plings have positive powers of w 2 /|q| 2 attached to them, 
with no other singular dependence on uo or q. Under 
z = 2 scaling, u scales to zero faster than |q|, and as a 
consequence all such couplings are irrelevant. 

Normally one would conclude that the critical behavior 
of the transition is mean-field like with logarithmic cor- 
rections, based on the fact that we are at the upper criti- 
cal dimension. While this may very well be the case here, 
it should be noted that the quadratic terms of the action 
(16) are quite unconventional, which may lead to uncon- 
ventional flow to the interaction u beyond tree level. We 
plan to study this as well as electronic properties of the 
system near the transition in future work. 

In an earlier work [11], Oganesyan, Fradkin and Kivcl- 
son obtained a bosonic theory for the transition from 
isotropic to nematic Fermi liquids by integrating out 
fcrmions via Hubbard-Stratanovich decoupling, similar 
to the Hertz-Millis theory. This transition corresponds 
to case m = 2 in our theory. These authors found z = 3 
and d c = 1. The origin of z = 3 lies in Landau damp- 
ing, which is a consequence of integrating out gapless 
fermions. In our approach there is no analog of Landau 
damping. Technically this is because we do not integrate 
out any gapless modes. More fundamentally however, it 
is clear that the transition is driven by a bosonic eigen 



mode whose energy goes through zero; by definition eigen 
modes are not damped. In the presence of interactions 
with other degrees of freedom of the systems, the wave 
function of this unstable mode needs to be determined 
self-consistently, and RG is a systematic procedure that 
allows the long wavelength mode to adjust itself to re- 
main an eigen mode as one approaches the low-energy 
limit. On the other hand if one integrates out the gap- 
less fermions up-front, no adjustment is allowed and the 
physically relevant modes appear damped. It thus ap- 
pears to be an artifact of this procedure. 
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